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Results are presented from high-precision computations of the orbital evolution and emitted gravita- 
tional waves for a stellar-mass object spiraling into a massive black hole in a slowly shrinking, circu- 
lar, equatorial orbit. The focus of these computations is inspiral near the innermost stable circular 
orbit (isco) — more particularly, on orbits for which the angular velocity fi is 0.03 ^ Sl/f2iaco < 1.0. 
The computations are based on the Teuksolsky-Sasaki-Nakamura formalism, and the results are 
tabulated in a set of functions that are of order unity and represent relativistic corrections to low- 
orbital-velocity formulas. These tables can form a foundation for future design studies for the LISA 
space-based gravitational-wave mission. A first survey of applications to LISA is presented: Signal 
to noise ratios S/N are computed and graphed as functions of the time-evolving gravitational-wave 
frequency for the lowest three harmonics of the orbital period, and for various representative values 
of the hole's mass M and spin a and the inspiraling object's mass fi, with the distance to Earth 
chosen to be To = 1 Gpc. These S/N^s show a very strong dependence on the black-hole spin, as 
well as on M and fi. Graphs are presented showing the range of the {M, a, fj,} parameter space, 
for which S/N > 10 at ro = 1 Gpc during the last year of inspiral. The hole's spin a has a factor 
~ 10 influence on the range of M (at fixed n) for which S/N > 10, and the presence or absence 
of a white-dwarf-binary background has a factor ~ 3 influence. A comparison with predicted event 
rates shows strong promise for detecting these waves, but not beyond about IGpc if the inspiraling 
object is a white dwarf or neutron star. This argues for a modest lowering of LISA's noise floor. A 
brief discussion is given of the prospects for extracting information from the observed waves. 



PACS numbers: 04.30.-hx, 04.80.-hz, 97.60.Lf 



I. INTRODUCTION AND SUMMARY 

Earth-based gravitational-wave detectors operate in 
the high-frequency band, ~ 1 — 10^ Hz, in which lie 
the waves from black holes of masses ~ 2 — 
Space-based detectors operate in the low-frequency band, 
~ 10~^-1 Hz populated by waves from black holes of 
mass ~ IO'^-IO^Mq. The high-frequency band is likely 
to be opened up early in the next decade by the LIGO- 
VIRGO network of earth-based detectors |^]. The pre- 
mier instrument for the low-frequency band is the Laser 
Interferometer Space Antenna (LISA) 

The European Space Agency has selected LISA as one 
of three "Cornerstone" missions in its "Horizon 2000-I-" 
program, NASA has appointed a mission definition team 
for LISA, and the ESA and NASA teams are talking to 
each other informally about the possibility of flying LISA 
as a joint ESA/NASA mission in the ^ 2010 time frame. 

One of the most interesting and promising gravita- 
tional wave sources for LISA is the final epoch of in- 
spiral of a compact, stellar-mass object into a massive 
black hole. In the LISA frequency band, where the 
central hole must have M ^ lO^M©, all giant stars 
and main-sequence stars will be tidally disrupted before 
the end of their inspiral, but compact objects — white 
dwarfs, neutron stars, and small black holes — can sur- 
vive intact. (Depending on the hole's spin, a massive 



white dwarf will be disrupted before the end of inspiral 
if M < Mmax ~ IC-IO^Mq. Neutron stars and small 
black holes can never be tidally disrupted in the LISA 
frequency band.) 

Sigurdsson and Rees |^ have estimated the event rate 
for such compact objects to spiral into massive black 
holes. "Assuming most spiral galaxies have a central 
black hole of modest mass (~ IO^Mq) and a cuspy 
spheroid," and for "very conservative estimates of the 
black hole masses and central galactic densities," they 
estimate one inspiral per year within IGpc distance of 
Earth. Most of the inspiraling objects are likely to 
be white dwarfs or neutron stars; the inspiral rate for 
stellar-mass black holes (/i ~ 6-10Mq) may be ten times 
smaller, about 3 per year out to 3Gpc, according to Sig- 
urdsson Q. Sigurdsson notes, however, that that the ev- 
idence for a recent burst of star formation in the central 
region of our galaxy suggests that normal nucleated spi- 
rals might have such starbursts every ~ 10^ years, which 
would enhance the stellar-mass black-hole density by a 
factor ~ 10 and would lead to stellar-mass black-hole in- 
spirals of one per year out to one Gpc. He notes, further, 
that if there was just one 50Mq black hole in the core 
of each galaxy now containing a ~ 10^ Mq central black 
hole, the result would be several inspirals of such 5OM0 
holes per year out to a cosmo logical redshift z = 1, all 
readily observable by LISA. 
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LISA'S observations of waves from such inspirals will 
have major scientific payoffs 

• Ryan |^ has shown that for circular equatorial orbits, 

the waves will carry, encoded in themselves, a map 
of the vacuum spacetime metric of the central black 
hole (or, equivalently, the values of the hole's multi- 
ple moments), and he has made a first, very crude, 
estimate of the precision with which LISA can ex- 
tract that map 0. Ryan's estimate is quite promis- 
ing. From the extracted map, one can determine 
whether the hole's geometry is that of the Kerr met- 
ric (i.e. "test the black- hole no hair theorem"), and 
one can use such maps to search for other kinds 
of conjectured massive central bodies (e.g. soliton 
stars il and naked singularities). It seems likely 
that this is true not only for circular geodesic or- 
bits, but also for generic orbits. 

• The observation of many such events will provide (i) 

a census of the masses and spins of the massive 
central holes, (ii) a census of the masses of the in- 
spiraling objects (which depend on and thus tell 
us about the initial stellar mass function and mass 
segregation in the central parsec of galactic nuclei), 
and (iii) a census of event rates (which depend on 
physical processes and on gravitational potentials 
in the central parsec). 

• In active galactic nuclei, the inspiral orbit may be sig- 

nificantly affected by drag in an accretion disk, 
producing both complications in the interpretation 
of the observations and opportunities for learning 
about the disks' mass distribution 

In planning for the LISA mission, it is important to 
understand the details of the waves emitted by such in- 
spirals. Those details are the most important factors in 
the choice of the mission's noise floor and its duration, 
and are likely to be the principal drivers of its data anal- 
ysis requirements and algorithms. 

The foundations for computing the emitted waves are 
nearly all in place: 

• If the orbit is known, then the waveforms and strengths 

can be computed using the Teukolsky - Sasaki- 
Nakamura ||ll| (TSN) formalism for first-order per- 
turbations of Kerr black holes. 

• The orbital evolution is governed by radiation re- 

action (and, if there is a robust accretion disk 
present, by accretion-disk drag j|]). Most massive 
holes are in galaxies with normal (non-active) nu- 
clei, and are thought to be surrounded by tenuous 
disks with "advection-dominated accretion flow" 
(ADAF). Narayan [|2j has shown that accretion 
drag should be totally negligible in such ADAF 
disks, so the orbital evolution is very cleanly gov- 
erned by radiation reaction. This is the situation 



that we analyze in this paper; we ignore accretion- 
disk drag. Those few holes that are in active galac- 
tic nuclei may be surrounded by "thin" or "slim" 
accretion disks, for which Chakrabarti and col- 
leagues Q have shown that accretion-disk drag may 
be significant. 

• The radiation reaction's influence on the orbit can be 

characterized fully by the rates of change of three 
"constants" of the orbital motion: the orbital en- 
ergy E, axial component of angular momentum 
L, and Carter constant Q |T^. From the emit- 
ted waves (computed via the TSN formalism) , one 
can read off _B = dE/dt and L; but the only known 
way to compute Q is directly from the radiation 
reaction force. 

• A formal expression for the radiation reaction force has 

been derived recently by Mino et. al. |l4j and by 
Quinn and Wald and several researchers are 
now working hard to convert this into a practical 
computational tool for deducing Q iQ. This will 
complete the necessary set of tools for computing 
all details of the emitted waves. 

The emitted waves will be so complex and so rich in 
structure and in parameter dependence, that it will re- 
quire extensive computations to give us the full knowl- 
edge required by the LISA mission. Those computations 
are proceeding in stages: 

1. Initial quick surveys, based on the Newtonian 
or quasi-Newtonian orbits and the quadrupole- 
moment approximation to gravitational-wave emis- 
sion. Such surveys are the foundation for the event 
rate estimates by Sigurdsson and Rees discussed 
above. 

2. More detailed and accurate surveys for orbits in 
the massive hole's equatorial plane, using the TSN 
formalism. Such surveys do not require computing 
Q, since Q vanishes for equatorial orbits. These 
surveys are of several types: 

(a) Studies of the evolution of the orbit's eccen- 
tricity. Such studies have been carried out by 
Tanaka et. al. and Cutler et. al. Il^] 
for non-spinning holes, and by Kennefick |2G] 
for small eccentricities around spinning holes. 
These studies, coupled with estimates of the 
orbital eccentricities when the objects are far 
from the hole and are being frequently per- 
turbed by near encounters with other objects 
piP], suggest that, despite the circularizing 
effect of radiation reaction, the eccentricities 
will still typically be large, e > 0.3, when the 
object nears the hole's horizon. 

(b) Systematic computations of the details of 
the emitted waves and the orbital evolution 
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for circular, equatorial orbits. This paper 
presents such computations and a compan- 
ion paper [ p2[ extends them to the transition 
regime, near the innermost stable circular or- 
bit (isco), during which the orbit makes a 
gradual transition from adiabatic inspiral to 
a plunge into the hole. 

(c) Computations of the waves' details and orbital 
evolution for elliptic, equatorial orbits. First 
explorations have been carried out by Shibata 
for general ellipticity and by Kennefick for 
small ellipticity (though in the 1970s and 
1980s there were studies for equatorial orbits 
that plunge from radial infinity into a hole or 
scatter off a hole . 

3. Surveys of the orbital evolution and waves for cir- 
cular orbits out of the hole's equatorial plane. It 
is known that radiation reaction drives circular 
orbits into circular orbits, thereby causing Q to 
evolve in a manner that is fully determined by TSN- 
formalism calculations of E and Q |25-2^]. There- 
fore, the tools are fully in hand for these surveys, 
and Hughes is in the late stages of the first 
one. (See Shibata ||2^ for an exploration of the 
wave emission before anyone knew how, correctly, 
to compute the orbital evolution, and see Shibata 
et. al. |3^] for studies of orbits with very small in- 
clination angles to the equatorial plane.) 

4. Surveys of orbital evolution and waves for the 
generic, most realistic situation: elliptic orbits out- 
side the equatorial plane. Such surveys must await 
a practical computation technique for Q. 

For circular, equatorial orbits (the subject of this pa- 
per), there have been extensive previous calculations, 
beginning with the pioneering study by Detweiler 
for a review see Mino et. al. jlj]. However, these pre- 
vious calculations have been motivated by the needs of 
LIGO/VIRGO observations in the high-frequency band, 
where (i) the ratio n/M of object mass to hole mass is 
not very small, so finite-mass-ratio effects (omitted by the 
TSN formalism) are important, and (ii) almost all of the 
observed inspiral signal comes from radii large compared 
to the hole's horizon, so post-Newtonian techniques can 
be used. The previous calculations have focused almost 
entirely on carrying the post-Newtonian calculations to 
very high order, on developing techniques for accelerat- 
ing their convergence, and — via comparison with TSN 
calculations — on evaluating their convergence jl^ . 

LISA'S regime and needs are quite different from this. 
For LISA, most of the signals are likely to come from 
systems with extreme mass ratios, /x/M <ti 1, for which 
(a) the TSN formalism is highly accurate and (b) the 
object lingers for a very long time in the vicinity of the 
hole's horizon before plunging into it. This means that 
post-Newtonian calculations are neither needed, nor ap- 
propriate. 



Because of these differences between the LIGO- VIRGO 
regime and the LISA regime, the previous TSN-based 
studies do not serve LISA's needs. The purpose of this 
paper is to begin fiUing that gap. Specifically: 

In this paper we introduce a new set of functions A/", 
T, 5, to characterize the orbital evolution and the emit- 
ted waves. These functions are dimensionless and of or- 
der unity, and depend on the hole's dimensionless spin 
parameter a = (angular momentum)/M^ and on the or- 
bit's dimensionless radius f = r/M. We give extensive 
tables of these functions, as computed by one of us (LSF) 
using the TSN formalism. We then use those tables to 
compute the evolution of the waves' frequency and sig- 
nal strength in LISA for a number of instructive values 
of the parameters M = (hole mass), a = (hole spin pa- 
rameter), /I = (object mass), and m = (wave harmonic 
order) = (wave frequency) /(orbital frequency). From 
these computations we draw a number of conclusions of 
importance for the LISA mission. 

The paper is organized as follows. Our notation, in- 
cluding the dimensionless functions Af, T, is intro- 
duced in Sec. Formulas for computing the dimen- 
sionless functions, and formulas for the orbital evolution 
and the waves' properties are given in Sec. IH. Tables 
of the dimensionless functions are given and discussed in 
Sec. Applications to LISA are presented in Sec. 0. 
Finally, concluding remarks are given in Sec. 0. 



II. NOTATION 

In this paper we shall adopt the following notation to 
describe the compact object's inspiral and the gravita- 
tional waves it emits; throughout we use geometrized 
units, i.e. we set G = (Newton's gravitation constant) — 
1 and c = (speed of light) = 1. 

fi: The mass of the inspiraling object. 

M: The black hole's mass. 

rj = fi/M: The mass ratio, assumed <C 1. 

a = S/A'P: The hole's "rotation parameter;" here S is 
the hole's spin angular momentum. 

r: The orbi t's Boyer-Lindquist radial coordinate; defined 
by ^r'^ + a'^{l + 2M/r) = (l/27r)x(the object's 
orbital circumference) 

tilde: A tilde over a quantity means that it has been 
made dimensionless by multiplying by the appro- 
priate power of M and, when the quantity is cx ^, 
multiplying by a factor l//i. 

r = r/M: The dimensionless radius of the orbit. 

il: The object's orbital angular velocity, as measured us- 
ing Boyer-Lindquist coordinate time t (defined be- 
low), i.e., using clocks that are far from the hole 
and at rest with respect to it. 
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= Mfl: The dimensionless orbital angular velocity, 
which is related to f by f2 = l/{f^^^ + a); cf. Eq. 
(2.16) of Ref. m. When Cl is small (large f), Ke- 
pler's laws dictate that ft ~ (M/r)^/^ = (orbital 
velocity)'^, i.e. (orbital velocity) ~ fl-^^^. 

subscript isco: A quantity evaluated at the object's in- 
nermost stable circular orbit ("isco"), where the 
inspiral ends and the plunge begins; for example, 
fiisco is the value of fl at the isco. 

t: Boyer-Lindquist coordinate time, or equivalently time 
as measured at "radial infinity" or on Earth. 

T: The Boyer-Lindquist time At until the isco is reached; 
i.e., the total remaining duration of the inspiral. 

-^orb^ The number of orbits remaining until the isco is 
reached. 

Tq: The distance from the binary to Earth. 

m: The order of a harmonic of the orbital frequency. 

fm = (rn/27r)f2: The frequency of gravitational waves in 
the m'th harmonic. 

E: The object's total energy including rest mass, i.e., 
the component —pt of its 4- momentum. Note that, 
because the object is gravitationally bound to the 
black hole, E < fj,, its gravitational binding energy 
is ^ - > 0. 

E = E/fi. 

Eoc- The total rate of emission of energy into gravita- 
tional waves that go to infinity. 

Eh- The total rate of emission of energy into gravita- 
tional waves that go down the horizon. 

Egw = Eoc + Eh ■ The total rate of emission of energy 
into gravitational waves that go both to infinity 
and down the hole's horizon; and also, by energy 
conservation, the rate of decrease of the object's 
total energy; i.e., E'gw = dEQYi /dt = —dE/dt. 

Eoom- The total rate of emission of energy into the m'th 
harmonic of the waves that go to infinity. 

ho.m = {hm+^ + h„ix^)'- The rms amplitude of the 
gravitational waves in harmonic m emitted toward 
infinity, at a time when the wave frequency is fm] 
here hm+{t,n) and hmx(t,n) are the two wave- 
forms emitted in a direction n and arriving at the 
Earth's distance Tq; (. . .) is an average over n and 
over a period of the waves; and the average over 
time automatically produces a factor 1/2 thereby 
making h"^ ^ be the mean value of ^ [(amplitude of 
hm+)^ + (amplitude of ft.m,x)^]- 



hc,m = ho.iny '^fm/ fm ■ A characteristic amplitude for 
the waves in harmonic m; here and throughout this 
paper the dot denotes a time derivative. The sig- 
nificance of hc^m is discussed below. 

K,m = /ic,m min[l, ^3(1 - /m//m,isco)]: A modified 
characteristic amplitude, discussed below. 



hn{f)= yfSf^^if) ■■ LISA'S "sky-averaged" rms noise 
in a bandwidth equal to frequency /. Here 
Sf^{f) is the one-sided spectral density Sh^{f) 
for some linear polarization 4-, inverse averaged 
over source directions and polarization ("sky- 
averaged"), l/S^^ = {1/3,,^). 

From the general relation di?Gw/'^^c?^ = (l/167r)(/i^-|- 
h^, ) for the energy flux in gravitational waves in terms 
of the time derivatives of the two waveforms (Eq. (10) 
of Ref. Il^), we infer that the rms amplitude and the 
energy in harmonic m are related to each other by 



ho,rn = . (2.1) 

We shall use, as our measure of where the object is 
in its orbit, the dimensionless orbital angular frequency 
f2, which is related to the gravitational-wave frequency 
in harmonic m by = (m/27rM)r2. We shall write 
various fully relativistic, time-evolving quantities (i?GW, 
hc.rm etc.) as the leading-order ("Newtonian") term in 
an expansion in fl^^^ ~ (orbital velocity), multiplied by 
relativistic corrections. Our notation for the relativistic 
corrections will be the following: 

Af: The correction to fP/Ct = CP/ft, where the dot is 
a time derivative. Note that ^P/Cl — d^/dlnfl is 
the number of radians d^ of orbital motion required 
to produce (due to radiation reaction) a fractional 
change dQ/Q in the orbital frequency. 

Afovh- The correction to A^orb (the number of orbits re- 
maining until the end of the inspiral). 

T: The correction to T (the remaining time to the end 
of the inspiral). 

£: The correction to -Eqw (the total energy loss rate). 

£oom'- The correction to Eoom (the energy radiated to 
infinity in harmonic m). 

Ho.m- The correction to ft-o,™ (the rms wave amplitude 
in harmonic m). 

'Hc.m' The correction to /ic.m (the characteristic ampli- 
tude in harmonic m). 
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The characteristic amplitude hc^m needs some expla- 
nation. As the object spirals inward in its orbit, its m'th 
harmonic waves spend ~ fm/fm = rf^m/(27r(iln/m) 
cycles in the vicinity of frequency (where $m is 
the harmonic's phase). Correspondingly, in a detector 
that observes the waves throughout the inspiral epoch 
Ai — fm/fm, the signal is enhanced, in comparison to 
the detector noise, by approximately the square root of 
this quantity. The signal strength is thus approximately 
the same as would be produced by a broad-band burst of 

amplitude hc,m = ho^m\J '^fml fm- 

The factor 2 inside the square root arises from a more 
precise definition of ft.c,m [0: The signal to noise ratio 
produced by the waves' m'th harmonic, averaged over all 
possible orientations of the source and the detector, is 
given by 



hc,m(^f 17 



1 2 



Kifm) 



d In /„ 



(2.2) 



where hn{fm) is the detector's rms noise at frequency /„, 
in a bandwidth equal to frequency, averaged over the sky. 
Equation (2.2) serves as a definition of /ic,m(/)- The re- 



lation hc^mifm) 
Eq. (29) ofRef. | 



{fm)^j2f^/ fm then follows from 
}] (with the factor 2 changed to 4 to cor- 



rect an error ), to gethe r w ith the definition of ho.m given 



above, Eqs. (2_^) and (2J^), and the evaluation of Fourier 
transforms using the stationary phase approximation. 

When the inspiraling object nears the isco, the band- 
width available for building up its signal in the de- 
tector becomes less than A/ = /. A good mea- 
sure of this reduced bandwidth is A/ = 2(/ — /isco) 
(with half of this band below / and half above). This 
is less than / for 2/isco/3 < / < /isco- Corre- 
spondingly, the amplitude of the built-up signal is ~ 
^0^2 (/isco — /)/(2/isco/3). Our modified characteristic 
amphtude h[ „^ = /ic,mniin[l, ^3(1 - /,„//,„, isco)] takes 
this signal reduction into account. 



III. FORMULAS FOR INSPIRAL AND WAVES 

In this section we shall give leading-order (in fl^/"^) for- 
mulas for the various time-evolving quantities, as func- 
tions of the dimensionless orbital angular frequency Ct 
and black-hole spin a, and thereby we shall produce ex- 
act definitions of the relativistic correction functions. To 
make clear the magnitudes of various quantities, we shall 
write some of our formulas numerically in a form rele- 
vant to LISA (for which we choose as a fiducial frequency 
/2 = O.OlHz and as a fiducial source, a = IOMq black 
hole spiraling into a M — lO^M© hole at ro = 1 Gpc 
distance from Earth). We shall also write our formu- 
las in a form relevant to the LIGO- VIRGO network of 
high-frequency detectors (with, as our fiducial frequency. 



lOOHz, and our fiducial source, a IMq neutron star spi- 
raling into a IOOM0 hole at IGpc distance). 

In the Newtonian limit, the orbital radius and orbital 
angular velocity are linked by the Keplerian relation 



r r 



(3.1) 



This permits us to write the number of orbital radians 
spent near orbital angular frequency in the following 
form (cf. Eqs. (3.16) of MTW H, in which a is our 
orbital radius r): 



n dinfi 9677175/3 

1.17 xlO^ /IOM0 
~ (/2/.01Hz)5/3 

_ 117 / IMq 

^ (/2/100Hz)5/3 [~jr 



M 
IOOM0 



M 



Af (3.2) 



Here Af is the general relativistic correction, which is 
unity in the "Newtonian" limit ^ 1. Similarly, 
the total remaining time until the end of the inspiral 
(Eq. (36.17b) of MTW) is 



5 1 M ^ 

T = ^ — r 

256 Tj 

1.41 X lO^^sec flOMQ\ flO^MQ 



(/2/.01Hz)8/3 \ ^ J \ M 
0.141sec /1Mq\ /IOOM0 



(/2/100Hz)8/3 \ n J \ M 



T 



T 



(3.3) 



and the number of orbits remaining until the end of the 
inspiral is 



iVorb = 7^ 
ZTT 



lnO_ d$ ^, ^ 111, 

a In S 2 = — _ ,„ A/orb 



InSl 



d\nn 



1.11 X 10"* 



lOMp 



647r f] 
M 



AAorb (3.4) 



(/2/.01Hz)5/3 ^ ^ 

" (/2/100Hz)5/3 [ )[ M ) '''"'^ ■ 

By integrating Eq. ( |3.2| ) inward to the isco, one can 
derive the following expression for the general relativistic 
corrections T for T and A/'orb for A^orb in terms of that 
M for VL^/h: 



3 Jn 



A/'orb = In'/' 



iiii/3' 



Mdn 



(3.5) 



The total energy loss rate [Eq. (3.16) of MTW] is 



(3.6) 



5 



5 



(3.7) 



where £ is the general relativistic correction. 

When the object is at large radii (small SI), the power 
radiated by the system's mass multipole moments /''^™ 
is of order ?7^S7^+^'/'^, while that radiated by its cm'rent 
multipole moments 5''=^™ is of order 77^ j^2+2(/+i)/3 ^ 

Correspondingly, the power Eoom radiated to infinity in 
harmonic m comes almost entirely from the moments of 
lowest allowed orders I, with the current moments of or- 
der I being, a priori, comparable to the mass moments of 
order / + 1. 

For TO = 1, the lowest allowed order for either mass 
or current is quadrupolar, since gravitational waves are 
always quadrupolar or higher. For circular orbits in the 
equatorial plane, the m = ±1 components of the mass 
quadrupole moment vanish, so the dominant waves are 
current quadrupolar 5'^''''^ and mass octupolar P-^^ . All 
other multipolar contributions to Eaoi are smaller than 
these by at least (flr)^ = fP^^. The contributions of 
5^'^^ to the radiated power Eooi can be derived from 
Eqs. (4.16), (5.27), (2.18a,c), (2.7) and (2.8) of Ref . |6|; 
and those of from Eqs. (4.16), (5.27), (2.7), and 

(2.8) of ll^. The sum of these dominant contributions is 



Eool = —if ^'^£001 
zo 



(3.8) 



where we have tacked on the general relativistic correc- 
tion factor £ooi to account for contributions from all the 
higher-order multipoles and to make the formula be valid 
not just for large orbital radii r but for all r > risco- In 
Eq. (U|), the numerical factor is 5/28 = 8/45 -t- 1/1260, 
where the big piece 8/45 is current quadrupolar, while 
the tiny piece 1/1260 is mass octupolar. 

For harmonic to > 2, the lowest allowed multipoles are 
of order I = to, and the mass moment /'"■i™ is nonzero 
so it dominates. All other multipolar contributions to 
Em are down from these by at least f]^/^ . An expression 
for Ern can be derived from Eqs. (4.16), (5.27), (2.7), and 
(2.8) of Ref. !§. The result is 



E„, — 



2(to + 1)(to + 2)(2to + 1)!to 



2m+l 



(to - 1)[2™to!(2to 

X„2Q2+2m/3c 



1) 



1112 



(3.9) 



where Soom is the relativistic correction and (2to + 1)!! = 
(2to-H1)(2to-1)(2to-3) • • -1. For to = 2, 3, and 4, this 
expression reduces to 



p _ 32„2ol0/3c 



E: 



- 243 2q4c; „ 
3 — 28 ' '-'oo3? 



(3.10) 



E, = ^^2j^l4/3^^^_ 



Note that the low- SI limit of E2 is identical to that of the 
total energy loss ii^cw [Eq. (|3.7D], as it must be since the 
TO = 2 harmonic dominates at low orbital velocities. 



From Eqs. (3.8) and ( |3.9[ ) for Em and the general rela- 
tionship (2T) between the waves' amplitude and energy, 
we obtain the following Newtonian-order expression for 
the amplitude in harmonic to: 




7 To 



(3.11a) 



m — 1 



TO -I- 1)(to + 2)(2to -I- l)!m2 
(to - 1)[2'"to!(2to-|- 1)!!]2 

^H^Qrn/Sj^^^ for TO > 2, (3.11b) 

ro 

where the relativistic correction is related to that for the 
energy by 



'^0,711 — '\J ^oo7n- (3.12) 

For the dominant, 7n = 2, radiation Eq. ( ^.llb|) becomes 

V 5 To 

3.6 X 10-22 f n \ f M Y ( f2 



rJlGpc \\QMqJ VlO^Afo 

3.6 X 10-23 / M W M 
r-o/lGpc \M^) V lOOM, 



.OlHz 



lOOHz 



(3.13) 



From Eqs. (|3.11 



for hn 



the definition of hr m m terms 



of ho,m, the relation fm — {m/2T:)fl, and Eq. (3.2), we 
obtain the following expression for the characteristic am- 
plitude in harmonic to 



hc.i 



V672^ 



(3.14a) 



/ 5(m +l){m + 2)(2m + l)\im?™- 
I 127r(TO- 1)[2"to!(2to + 1)!!]2 



fj(2m-5)/6-^^^^ for m > 2, (3.14b) 



where [using Eq. (3.12)] the relativistic correction is re- 
lated to earlier ones by 



He 



(3.15) 



For TO = 2, expression ( 3.14b ) becomes 
"2"r?i/2M, 



hn2 = 



37r roS]i/6 
1.0 X 10-1^ 



ro/lGpc \IQMqJ \1Q^Mq 



( M 



3.2 X 10-22 / ^ 
ro/lGpc \M^^ 



.OlHzV 



M y /lOOHzV 



c,2 



(3.16) 
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All of the relativistic correction functions can be ex- 
pressed analytically in terms of £, and foom- This has 
almost been done already: The correction functions T, 
'Ho,m and Hc.m have been expressed in terms of £, Soom, 

and tlisco by Eqs. (3^), ( 3.12| ) and ( 3.15 ) respectively. 
All that remains is to derive an expression for M in terms 
of £, and an expression for ilisco- 

The derivations are based on the Kerr-metric relations 



E : 



1 - 2/?^ + a/r3/2 
Vl-3/f + 2^V^ 



(3.17) 



for the object's total energy in terms of its dimensionless 
orbital radius f (Eq. (5.4.7b) of Ref. E^), and 



f = (Q-^ - a) 



2/3 



(3.18) 



for its orbital radius in terms of its orbital angular veloc- 
ity (Eq (2.16) of Ref. [|2|). By differentiating these equa- 
tions with respect to tim e an d combining with each other 
and with Eqs. (HJ) and (jsj) for Q'^/h and i;Gw = -E, 
we obtain 



£ 



a 


y/3 


f3/2 




3 


2a 


- + 
r 


^3/2 



1-^ 

r 

-3/2 



8a 

f3/2 



3a^ 



(3.19) 



When f is regarded as the function ( |3.18 ) of tl, this be- 



comes the desired expression for Af in terms of £ and 

n. 

The innermost stable circular orbit (isco) is at the lo- 
cation fisco where the object's total energy E(f) is a min- 
imum, or equivalently where O is infinite, or equivalently 
where Af vanishes; i.e., fisco is that root of the quartic 
equation — 6f + 8af — 3a^ = which lies between 
1 (when a = 1) and 6 (when a = 0). An analytic ex- 
pression for Tisco has been given by Bardeen, Press, and 
Teukolsky M: 



fisco = 3 + Z2 - sign(a)[(3 - Zi)(3 + Zi + 2Z2 



1I/2 



l + (l-a2)i/3[(i + a)i/3 + (i_„)i/3] ^ 
{3a^ + Zi^y^^. (3.20) 
The dimensionless orbital angular velocity at the isco 



Zi 
Z2 



ilisco is expressed in terms of this fisco by Eq. (3.18) 



^isco 



1 



-3/2 



(3.21) 



We note, in passing, approximate analytic formulae 
for the relativistic corrections T and A/'orb to the time 
T and number of remaining orbits A/'orb until the end 
of inspiral — for the special case of a no nspinning black 
hole, a = 0. Inserting expression (|3.19|) into Eqs. (|3. 



and (3_^), and noting from Table ^ that in each of these 
equations £ is a much more slowly varying function of 



the integration variable than the rest of the integrand, 
we pull £ out of the integral (i.e., we perform the first 
step of an integration by parts) and then perform the 
integration analytically. The results are 



£ 



1 



2205, ,3 
T6~" ■ 



19845 ,, 4 
16 " 



A/'orb 



4671 
" £ 



19845 



3u 



32 



uMn 



3(l + \/2)2 



(1 + VT^^y 



Au 



288u3 



VI -3u 



(3.22) 



(3.23) 



where u = tl^/^ ~ 1/f. These formulae agree with the 
numerical values of T and A/'orb in Tables IX and ^ below 
to within 3 per cent at 6.02 < f < 18, and to within 1 
per cent at 6 < f < 6.02 and f > 18. 



IV. TABLES OF RELATIVISTIC CORRECTION 
FUNCTIONS 

We shall use two dimensionless parameters to measure 
the distance of an orbit from the isco: the ratio r/nsco = 
f/f isco of the orbit's Boyer-Lindquist radial coordinate r 
to its value at the isco, and the ratio rj/fJisco = r2/r2isco 
of the orbit's angular velocity to that at the isco. The 
relati onshi p betw een th ese t wo parameters is given by 
Eqs. (|3.18D , (|3.20|) and ( |3.2l| ), and is tabulated in Table 



We have integrated the Teukolsky-Sasaki-Nakamura 
equation for perturbations of a Kerr black hole to ob- 
tain the functions £(17) and £oom(^^), and we have then 
used Eqs. ( p^ , ( pT9| ), ( p^ ) and ^) to compute TV, 
T and A/'orb- These functions are listed in Tables llj-g 
and some of our numerical methods are described in the 
Appendix. As a byproduct of these calculations, we have 
inferred what fraction E/ Eq^ of the total rate of energy 
emission goes down the hole's horizon; that fraction is 
shown in Table VIE 



V. APPLICATIONS TO LISA 

A. LISA Noise 

Tentative error budgets for LISA are spelled out in 
Tables 4.1 and 4.2 of the LISA Pre-Phase-A Report 
Various researchers have computed LISA noise spectra 
from those error budgets. It is conventional, for LISA, 
to characterize the noise by the sensitivity to periodic 
sources for one year integration time and a signal-to- 
noise ratio of 5, averaged over source directions and po- 
larizations ( "sky averaged" ) . We shall denote this quan- 
tity by /igN5 lyr- ^ is related to the sky- averaged spec- 



tral density introduced in the paragraph before Eq. (2.1) 
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by /isN5,iyr = 5^^SA^^^ ^here A/ = 1/lyr is the 
bandwidth for the one-year integration time; and cor- 
respondingly, it is related to the sky-averaged rms noise 

in a bandwidth equal to frequency, hn{f) — \J fS^^{f) 
(which we use in this paper), by 



least To = IGpc. Accordingly, we have computed the 
range of masses fi and M and black-hole spins a for which 
(5/iV)ims > 10 at a distance ro — 1 Gpc. This range of 
"detectable systems" is shown in Fig. for LISA without 
the white-dwarf background (solid curves) and with the 
background (dashed curves). [Q 




.(/) • 



(5.1) 



We have deduced /i„(/) using this equation and the 
values of /i|N5,iyr(/) computed by various researchers 
p9|-pt; we plot it as a thick solid curve in Figs. |^-|^ 
below. 

It is likely that LISA's performance will be compro- 
mised at / < 0.003 Hz by a stochastic background due 
to white-dwarf binaries. The most recent estimate of 
that stochastic background is by Hils and Bender p2[ ; 
it agrees satisfactorily with an estimate by Webbink and 
Han We have used a simple piece- wise straight- 

line fit to the logarithm of the Hils-Bender white-dwarf- 
background noise curve: straight lines that join the fol- 
lowing points in (log^g /, log;^ 
sured in Hz: 



no 'ilNs.iyr) where / is mea- 



(-4, -20.518) , (-3.62, -20.737) , (-2.78, -21.66) , 
(-2.61,-22.90), (-2,-23.731). (5.2) 

This wh ite-dwarf noise, converted to our conventions via 
Eq. (5T), is shown as a thick dashed curve in Figs. 
below. 



B. Detectable Systems 

Of greatest interest, for probing the spacetime geome- 
tries of massive black holes, is the gravitational radia- 
tion emitted during the last year of inspiral of a compact 
object. In planning the LISA mission, it is important 
to know the detectability of these final-year waves, as 
a function of the system's parameters: the hole's mass 
M and spin a, the object's mass /i, and the distance ro 
from Earth. Previous studies of this issue [^,^ have 
assumed that the massive hole is nonspinning, a — 0. 

It is straightforward to compute the rms signal to noise 
ratio (5'/iV)rnis (aver age d over detector and system ori- 
entations) from Eq. ( |2.2| ), using the noise amplitudes hn 
described above and the dominant m — 2 characteris- 



tic amplitude he. 2 of Eqs. (p. 16) and ( 3.15 ), with M and 



^oom taken from Tables fVIII and IV. In this calculation 



the frequency /2 — I ^/tt ranges from its value at time 
T = I year [Eq. (Q and Table Q to its value at the 
isco. 

In view of the complexity of the data analysis for these 
waves, a signal to noise ratio of about 10 may be required 
for their detection, and in view of the estimated event 
rates (Sec. |), it is necessary that LISA see out to at 




FIG. 1. The minimum mass ^min that the inspiraling 
object must have in order to produce a signal to noise ratio 
S/N > 10 in its dominant harmonic, m = 2, during the last 
year of its inspiral. This /imin is plotted as a function of the 
black-hole mass M, for various black-hole spin parameters a. 
The solid curves are for the LISA noise spectrum; the dashed 
curves are for the LISA noise plus a stochastic background 
noise due to white-dwarf binaries, 

Several features of this figure deserve comment: 

• Inspiraling white dwarfs and neutron stars (/i ^ 
IAMq) are barely detectable, with (5/iV)i.ms = 10, 
at IGpc. It would be highly desirable to reduce 
LISA'S design noise floor by a factor two or three, 
to give greater confidence of detection. 

• For ji = lOM© inspiraling black holes, the de- 
tectable systems have a wide range of central black- 
hole masses: 10"* M„ < M < l{f Mo 



IQ. 



• The upper limit on detectable central-hole masses 
M depends strongly on the black-hole spin: for 
fj, = IOMq it ranges from 2 x 10*^ M© to 3 x IO^Mq 
without the white-dwarf background, and 1 x 10® 
to 1.5 x 10^ with the background. (The spins shown 
are for no rotation a = 0, and for near the maxi- 
mum rotation, a ~ ±0.998, that can be produced 
by spinup via accretion from a disk |^5[| .) 

• The white dwarf background reduces the maximum 
detectable black hole mass by about a factor 2.5, 
independent of the spin. 

• The white dwarf background and the black hole 
spin have little influence on the minimum de- 
tectable mass M. This is because, at low M the 
object travels a large radial distance in its last year 
of life, so most of the signal to noise comes from 
radii r 3> ^isco where the spin is unimportant, and 
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(by virtue of the small M) most comes from fre- 
quencies high enough that the white-dwarf back- 
ground is negligible. 

For probing the immediate vicinity of the hori- 
zon, we are interested in waves with frequencies, say, 
(2/3)/2,isco < /2 < /2,isco- Figure | shows the range of 
systems for which {S/N)^^s > 10 in this frequency band, 
at a distance ro = IGpc, during the last year of inspi- 
ral. Note that restricting attention to this near-horizon 
frequency range has reduced substantially the set of de- 
tectable systems: for fj, = 10 Mq, the minimum black-hole 
mass is increased by a factor ~ 20 to 100, depending on 
the spin a. Nevertheless, there is still a wide range of 
systems accessible for study. 




FIG. 2. The minimum mass /imin that the inspiraling 
object must have in order to produce a signal to noise ra- 
tio S/N > 10 in its dominant harmonic, m = 2, during the 
last year of its inspiral, and in the vicinity of the horizon, 
(2/3)/2,isco < fi < /2,isco. This ^min is plotted as a function 
of the black-hole mass M, for various black-hole spin param- 
eters a. The solid curves are for the LISA noise spectrum; 
the dashed curves are for the LISA noise plus a stochastic 
background noise due to white-dwarf binaries, 

For distances larger than r-Q ~ 1 Gpc, cosmological ef- 
fects have a significant influence on the signal |^6|. At 
fixed {a,M(l -I- z),iJ,{\ + z)} (where z is the cosmologi- 
cal redshift), the characteristic amplitude and signal to 
noise ratio scale oc I/tol, where Tol is the luminosity dis- 
tance to Earth. The scaling of (S'/A^)rms with n{l + z) is 
not so simple, because it influences the waves' frequency 
evolution in complicated ways that entail the relativistic 
correction functions. For extremely rough estimates, one 
can use the leading-order (in O) expression for hc.2 [Eq. 
( P6| ), S/N IX. hc,2 oc fJ,^''^/roh] to infer (l-l-z) cx rl^ 
for the minimum detectable object mass at flxed a and 
Af(l -I- z), but for reliable results, one must repeat the 
analysis (sketched above) by which we arrived at Figs. |l| 
and |. 



C. Evolution of the Waves During Inspiral 

To gain insight into the emitted waves and how they 
evolve during the inspiral, we have constructed Figs. |^— 
g. Each fi gure depicts the waves' evolution for the value 
of object mass n and hole mass M (in solar masses) listed 
in bold letters in the upper right corner. The horizontal 
axis is the waves' frequency / and the vertical axis, their 
modified characteristic amplitude h'^. As the inspiral pro- 
ceeds, the waves sweep upward in frequency (left to right) 
along one of the thin curves. These evolutionary curves 
are shown for three different values of the black-hole spin, 
a = —0.99 (retrograde orbit; short-dashed curves), a = 
(no rotation; long-dashed curves) and a = -1-0.999 (pro- 
grade orbit; solid curves). For each spin, three curves 
are shown corresponding to the three lowest harmonics 
TO = 1, 2, 3 of the orbital frequency. The values of a and 
TO for each evolutionary curve are listed near the vertical 
endpoint of the curve. Also shown in each figure is the 
rms noise amplitude /i„ for LISA: a thick solid curve in 
the absence of a white-dwarf-binary background and a 
thick dashed curve including that background. 

The range of frequency sweep is strongly dependent on 
the masses /i and M of object and hole. Neglecting the 
relativistic correction factor T (which is unimportant for 
this purpose when the frequency sweep is substantial), 

Eq. (P) tells us that /isco//lyr oc (^/M)3/8(l/M)3/8, 



where /lyr is the frequency one year before reaching the 
isco. Thus, the greatest frequency sweep is for the least 
extreme mass ratio and the smallest hole mass, fi/M = 
10/10^ (Fig. I) with /a sweeping from - 0.006 Hz to 0.4 
Hz; while the smallest sweep is for the most extreme mass 
ratio and largest hole mass, /i/M = 10/10^ (Fig- 0) with 
/2 sweeping only from - 0.0023 to 0.0027 Hz. 

The height of a signal curve h'^ above the noise curve 
hn is about equal to the signal to noise ratio in an ap- 
propriate bandwidth A/: A/ ~ f well away from the 
endpoint of inspiral, and A/ — 2(/ — /isco) near the end- 
point; cf. the discussion of the definition of h'^ at the end 
of Sec. |l| Near the endpoint of inspiral h'^ plunges for 
three reasons: (i) because of the narrowing of our cho- 
sen bandwidth; (ii) because the rate of frequency sweep 
speeds up due to flattening of the effective potential for 
the object's radial motion, and this produces a reduction 
in the number of cycles A'cyc in a given bandwidth and 
reduction in h'^ oc ^JN^; (iii) because, for large a and 
prograde orbits, the orbit sinks deep into the throat of 
the hole's embedding diagram, from where waves have 
difficulty escaping. 

On each signal curve there are three solid dots. They 
label (/, h'^ for specific times during the inspiral: T = 1 
year before the endpoint (leftmost dot), T = 1 month 
before the end (center dot), and T = 1 day before the end 
(right dot). Beside the dots for the dominant harmonic, 
TO = 2, are shown two numbers that characterize the 
orbit and waves at that time: the radius r — r / M oi the 
orbit in units of the black-hole mass, and the number of 
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gravitational- wave cycles in the m = 2 harmonic, from 
that time until the endpoint of inspiral. At the bottom 
end of each m — 2 curve is shown the radius f of the isco. 



o 




0.002 



1 CO — ^1 » 

Pis. 97 'lyp, "i° 



0.005 



0.01 
/ Hz 

FIG. 3. Gravitational waves from a IMq white dwarf or 
neutron star spiraling into a IO^A/q black hole at IGpc dis- 
tance from Earth, as observed by LISA. The thick solid curve 
is LISA'S rms noise level hn{f) averaged over the sky; the 
thick dashed curve is an estimate of the stochastic-background 
"noise" produced by white dwarf binaries. Each thin curve 
is the modified characteristic amplitude h'c{f) for a harmonic 
of the waves, and is labeled vertically by the hole's spin pa- 
rameter a and the harmonic number m. The three dots on 
each curve indicate the waves properties one year (left dot), 
one month (center dot) and one day (right dot) before reach- 
ing the isco. The dots on the dominant, m — 2, harmonics 
are labeled by the orbital radius r = r/M and the number of 
m — 2 wave cycles remaining until the isco. The isco radius 
is shown at the bottom of each m = 2 curve. 




0.001 0.002 



FIG. 4. Gravitational waves from a IOM0 black hole spi- 
raling into a 10^ M0 black hole at IGpc distance from Earth, 
as observed by LISA. For notation see caption of Fig. ^. 

It is worthwhile to scrutinize the details of these fig- 
ures, including the numbers beside the dots. Consider, 
for example, Fig. ^ for a fi = IOMq object (black 
hole) spiraling into a M = lO^M© hole. If the big 
hole is rapidly rotating and the orbit is prograde so 
a — -1-0.999, then the dominant m = 2 evolutionary curve 
shows the object, one year before its death, at f 6.80 



(3.4 Schwarzschild radii), with a signal to noise ratio of 
hc/hn ~ 100, and with 185,000 cycles of gravitational 
waves left until death. One month before death, the 
object is at f = 3.05 (1.53 Schwarzschild radii), with 
hc/hn ~ 50, and with 40,000 cycles left. One day be- 
fore death, it is at f = 1.30 (compared to 1.18 for the 
isco), with hc/hn ~ 10 and with 2,320 cycles left. It is 
impressive how long the object lingers in the vicinity of 
the horizon, and how many wave cycles it emits. 

For a nonspinning hole a = 0, the numbers are less 
impressive but still remarkable: the last year is spent 
spiraling from r — 9.46 (4.73 Schwarzschild radii) to the 
isco at f = 6 (3 Schwarzschild radii), during which 85,000 
wave cycles are emitted and hc/hn drops from ~ 100 to 
^ 10 at one day and then to zero. 

The large number of wave cycles carry a large amount 
of information about the source. We shall discuss this 



issue in Sec. VD below 
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FIG. 5. Gravitational waves from a IOOA/q black hole spi- 
raling into a 10^ M0 black hole at IGpc distance from Earth, 
as observed by LISA. For notation see caption of Fig. ^ 



Figures illustrate the influence of the mass of the 
inspiraling object on the signal strength. For M fixed 
at 1O^M0, one year before merger the m — 2 signal to 
noise ratios hc/hn are ~ 15 for /i = IM©, ~ 100 for /i — 
IOMq and 500 for = IOOMq. This is a moderately 



faster growth than our crude estimate oc p}!"^ in Sec. VB. 
Notice that hc/hn drops below 10 one month before the 
endpoint for /i = IMq, and one day before the end for 

/I = IOMq. 

To maximize the exploration of the horizon's vicinity, 
we want the object to spend its entire last year at radii 
f ^ 10. If the object is a IOA^q hole, this is the case when 
M > IO^Mq; cf. Figs. |,|6|and|. For M < IO^Mq, such 
exploration is debilitated by the large frequency sweep; 
cf. Fig. 1^. We have previously met this issue in Sec. VB. 

Figure ^ shows that the white-dwarf-binary back- 
ground is a serious issue for hole masses M ^ IO^Mq, 
while Figs. Ij-^show that it is relatively unimportant for 
M < IO^A/q. We have previously met this in Sec. VB 
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FIG. 6. Gravitational waves from a IOMq black hole spi- 
raling into a 10^ Mq black hole at IGpc distance from Earth, 
as observed by LISA. For notation see caption of Fig. H. 
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FIG. 7. Gravitational waves from a IOMq black hole spi- 
raling into a IO^A/q black hole at IGpc distance from Earth, 
as observed by LISA. For notation see caption of Fig. 0. 



D. Information Carried by the Waves 

As is well known the waves' highest accuracy in- 
formation is carried by the time evolution of their phase. 
For circular, equatorial orbits, where there is no orbital 
precession, the phase evolution will be the same for all 
the harmonics as for the orbit itself, and that phase evo- 
lution is embodied in d^/d\nil = 17^/0. Equation (3.2) 
shows this quantity, at fixed frequency, to be proportional 
where 



to AA/Ml,p, 



(5.3) 



is the system's chirp mass. Since a year of observations 
will typically entail A^cycic ~ 10^ cycles of waves, and by 
the method of matched filters one can detect a secular 
shift of one waveform with respect to another by a small 
fraction of a cycle [M , the "raw" precision for measuring 
the evolution of M/A'Ichivp will be of order 10^^. 

If most of the last year is spent near the horizon, say 
at frequencies ///isco — f^/^^isco ^0.1 (as will usually be 
the case), then this phase evolution will depend strongly 



not only on the chirp mass, but also — through the func- 
tion A/'(///isco) — on the black-hole spin parameter a. 




fif isco 

FIG. 8. Af the relativistic correction to d^/dlnO, = fl^ /Cl 
(the number of radians of orbital inspiral per unit logarith- 
mic change of orbital or gravitational- wave frequency), plot- 
ted against ///isco = n/fiisco (the ratio of gravitational- wave 
frequency to the frequency when the isco is reached and the 
inspiral ends). 

This strong a-dependence is exhibited in Fig. ^. Even 
for a < 0.5, where the curves A/'(///isco) for different a 
look very close together, dN /da ~ 0.1, this a-dependence 
translates into dNcydcs/ da ~ 10^, which is huge. Thus, 
it is reasonable to expect the measured phasing to de- 
termine both a and Mchirp to high precision — though a 
detailed parameter study is needed to be absolutely cer- 
tain. 

The absolute frequencies associated with the observed 
phase evolution (e.g., the measured frequency at the end 
of inspiral) are determined by a combination of a and the 
hole's mass M. This absolute frequency scale presumably 
will be measured much less accurately than the phasing 
itself, but still, probably, accurately enough to determine 
the mass Af to a very interesting precision. Knowing 
-M^chirp, a, and M, one can then compute the object's 
mass fi; and from the absolute amplitudes of the waves 
one can then infer the distance Tq from the system to 
Earth. 

Poisson has estimated the accuracies with which 
such phase-evolution measurements can determine A/, 
rj = fi/M, and a. His estimates are based on an ana- 
lytic model of the signal in which (translated into our 
notation) M is expanded in powers of Q — r^isco and 
only the leading order term is kept. Poisson assumes 
M = lO'^Afo, = 1OA//0, a ^ (i.e., not close to 
±1), and a measurement time of one year. For these 
parameters, our Figs. || and || suggest that, for a < 0.5, 
his expansion may be accurate to within a few tens of 
per cent, while for a ^ 0.9 it is seriously wrong. His 
estimated measurement accuracies are Aa ^ 0.05/p, 
AM/M - 0.002/p, and Srj/r] - 0.06/p, where p is the 
amplitude sig nal to noise ratio. 

Our Tables |VIIl| and |^ for TV and T (the relativistic 
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corrections to the orbital phase evolution rate d$ /d In / 
and the time T to the end of inspiral) can serve as the 
foundation a more definitive computation of the phasing- 
based measurement accuracies. 



0.14 




a=-.99 












0.12 














d 








=5_ 0.1 




//"//- y 




^ 0.08 




2.// ' , 




0.06 






7-^-^.9 
















0.04 
















0.02 






^^3e5:-^9 








^^~%g99 




0.1 


0.2 0.3 0.4 


0.5 0.6 0.7 



FIG. 9. The ratios hc,i/ha,2 and /ic,3/ftc,2 as functions of 
the black hole spin a and the orbiting object's radius r/risco- 

Information is also carried by the relative amplitudes 
of the waves' harmonics. Most promising, we think, are 
the amplitude ratios for the first and second harmonics 
and for the third and second. We plot these ratios in 
Fig. ^, as parametric functions of the hole's spin a and 
the orbital radius r/rjsco- From this plot it is evident 
that the instantaneous amplitude ratios will give both a 
and the instantaneous r/r^sco with moderate accuracy — 
though only for those systems with strong enough signals 
that the weakest of these harmonics, to = 1, stands up 
strongly above the noise; cf. the short-dashed curves in 
Figs.|^^. 

In our idealized case of circular, equatorial orbits, this 
harmonic-ratio information is not independent of that 
from the orbital phasing, but it could provide a confir- 
mation of the phasing conclusions. 

In the more realistic case of noncircular, nonequatorial 
orbits, the waveforms will be much richer and there will 
be many more parameters to solve for. Our survey of 
the circular, equatorial case gives some rough indication 
of the kinds of information one can extract and by what 
methods. 



VI. CONCLUDING REMARKS 

In this paper we have tabulated the results of TSN- 
based computations of the waves emitted by an object 
spiraling into a spinning, massive black hole on a slowly 
shrinking, circular, equatorial orbit. Our Tables |l|-^ 
can serve as a foundation for future mission-definition 
studies for LISA — most particularly, for studies of how 
changes in the mission design may affect LISA's ability 
to detect such inspiral waves, for studies of the accura- 
cies with which LISA's data can extract the properties of 
the source, and for explorations of possible data analysis 



algorithms. 

Much more important, in the long run, will be the ex- 
tension of our analysis to nonequatorial and noncircular 
orbits. This extension is urgent, since models of active 
galactic nuclei predict, rather firmly, that the orbits will 
be nonequatorial and quite noncircular, and since the ear- 
liest possible date for LISA to fly is less than ten years 
in the future. 
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APPENDIX: NUMERICAL METHODS 

Teukolsky [|o| found that the equations describing per- 
turbations of the Kerr spacetime could be separated into 
separate radial and angular equations. For the circular, 
equatorial orbits studied in this paper, the challenges of 
solving the perturbation equations are all associated with 
the numerical solution of the radial equation. We have 
used Green function methods to solve the radial equa- 
tion and determine the power radiated down the horizon 
and to infinity by a particle in a circular equatorial orbit. 
The general method of solution and formulation of the 
problem is well described in Ref. |49|, and we refer the 
interested reader there for details. In this appendix we 
describe several innovations that can dramatically speed 
the solution of the radial equation compared to the more 
conventional methods applied elsewhere. 

The Teukolsky radial equation is a second order, ordi- 
nary differential equation. In the form given originally 
by Teukolsky the equation is stiff and the solution 
satisfying the physical boundary conditions is difficult to 
obtain. Sasaki and Nakamura found, through a local 
change of variables, a form of the radial equation which 
is not stiff, and we have worked with the radial equation 
in that form. 

In the Sasaki-Nakamura formulation, the homogeneous 
(source-free) radial equation takes the form 



d , ^ d , ^ 
— -T{r,)—-U{r,) 
dr^: dvi. 



(Al) 
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Here is the so-called tortoise coordinate, which ranges 
from — oo at the (outer) horizon to oo at spatial infinity. 
(In this appendix, and in this appendix only, we express 
all dimensioned quantities, such as r and r*, in terms of 
the black hole's mass A/, eschewing for convenience the 
superscript-tilde notation used elsewhere in this paper.) 
The tortoise coordinate can be expressed analytically in 
terms of the Boyer-Lindquist radial coordinate r and the 
location of the inner and outer horizons r_|_ and r_ : 

2 

r* = r H [r_|_ In (r — r+) — r_ In (r — r_)] (A2a) 



where 



(r — r^){r — r^) — — 2r + 



and 



(A2b) 
(A2c) 



The functions J-' and U are parameterized by the angular 
frequency of the perturbation lu — 27r/, the angular mo- 
mentum of the spacetime a, and the angular separation 
constants £ and m (with |m| < i). (For the particular 
forms of T and U see ^^.) For circular, equatorial or- 
bits LU is always an integer multiple of the orbital angular 
frequency, uj = ujm = mfl = 27r/m . 

To obtain the Green function solution to the radial 
equation with source we need the two solutions to the 
homogeneous equation corresponding to the physical 
boundary conditions at infinity (no in-coming radiation) 
and the horizon (no up-going radiation). These solutions 
are determined numerically by posing the boundary con- 
ditions near infinity or the horizon and integrating the 
radial equation inward or outward, as appropriate. In 
the Sasaki-Nakamura variables, obtaining a solution to 
the radial equation poses no particular challenge; corre- 
spondingly, it is conventional to use a "work-horse" inte- 
grator {e.g., Runge-Kutta or Bulirsch-Stoer) to solve the 
equation. On the other hand, the radial equation arises 
from a separation of variables and is parameterized by 
the separation constants £, m and w, corresponding to the 
resolved angular and temporal dependence of the pertur- 
bation. Consequently, it is necessary to solve the radial 
equation separately for every important set of angular 
multipoles {£, m) and frequency uj. For very relativis- 
tic orbits even moderate accuracy in the total radiated 
power may require solving the radial equation tens of 
thousands of times for different angular multipoles and 
harmonics of the orbital frequency. Consequently, speed- 
ing the solution while preserving its accuracy is of fun- 
damental importance. In the remainder of this appendix 
we address several innovations we have made in solving 
this equation that, depending on the details of the orbit 
and the desired accuracy of the solution, can result in a 
several order of magnitude reduction in the solution time 
compared to a conventional approach. 



1. Boundary conditions at the horizon 

As one approaches the horizon, the physical solution 
for the radial function, corresponding to down-going ra- 
diation, leads to the boundary conditions used for the 
numerical integration of one of the homogeneous solu- 
tions of the Sasaki-Nakamura equation: 



lim Xnir^,) = e 

— > — oo 

hm — — (r* ) = - 



where 



UJ- 



am 
2^ 



(A3a) 
(A3b) 

(A3c) 



and r+ is the radius of the outer horizon in Boyer- 
Lindquist coordinates. 

As a practical matter the boundary conditions used to 
determine Xh are posed at some large, negative but finite 
r*, sa y ; i.e., "close to", but not at, the horizon. Using 
Eqs. (A3) evaluated at finite i?l for the boundary condi- 
tions introduces errors of fractional order (S = r — r+ into 
the solution. This error can be represented as an error 
in the amplitude of the power radiated down the horizon 
and the introduction of some small component of radia- 
tion up-going from the horizon. These errors propagate 
to large where they contribute to the out-going radia- 
tion and lead to errors in the calculated power radiated 
to infinity by the orbiting particle. 



The errors introduced by using equations |A3| when pos- 
ing boundary conditions at finite radius can be expressed 
as a power series in S. The coefficients of that expansion 
can be estimated by solving the equations several times, 
for different R'H , and using Richardson extrapolation. To 
estimate the first N terms in the error expansion requires 
iV + 1 numerical solutions of the equations, each begin- 
ning with the boundary conditions posed at a different 
R'^. Controlling the error requires that the radial equa- 
tion be solved at least twice and often three or more times 
at different, large 

To improve the convergence rate of this error estimate 
and allow us to pose our boundary conditions at smaller 
we have solved the Sasaki-Nakamura equation an- 
alytically about the point at = — oo, finding the first 
corrections in S to the boundary equations given by Eqs. 
([A3|) . The improved boundary conditions are given by 



X{n) = 1 + 5A' 



5d 



-iuj^ (1 + 5A') + A' 

2r+ 



(A4a) 
(A4b) 



where 



A' 



^(1) 



4ri 



(1) 



G' 



-f it LO- 



(A4c) 
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S : 

d. 



U. 



(1) _ 



+ J 

d 
2^ 



4r_ 

A H 2r+a)_ 



2(2 + d) 2r^ 



d2 



- 1- 



rf(2 + d)\ 4r 



d^ 



(A4d) 
(A4e) 

(A4f) 



(A4g) 

(A4h) 
(A4i) 
(A4j) 

(A4k) 

(A41) 
(A4m) 

(A4n) 

(A4o) 
(A4p) 
(A4q) 



The numerical solution to the radial equation using 
these improved boundary conditions converges upon the 
true solution more quick ly than a solution using the 
boundary conditions (A3). We are thus able to pose ap- 
proximate horizon boundary conditions at smaller \R*L\, 
reducing the domain over which we must integrate the 
radial equation and, often the number of times we must 
integrate the equation for each (Lu,i,m) in order to ob- 
tain a solution of controlled accuracy. 



2. Boundary conditions at spatial infinity 

As r« ^ oo, the physical solution for the radial func- 
tion, corresponding to no in-going radiation, leads to the 
boundary conditions for the numerical integration of the 
other critical solution of the radial equation: 





= 2ir 


+c^_F(0) , 




= -2 


+ d/2. 


F(0) 




7(r) 








k=() 




70 


= A (A -)- 2) — 12aa; {au — m) — \2iijj , 


71 


EE 8i [ 


Sa^u! — aX {aoj — m)] , 


72 


= 12 


—2ai (au! — m) + — 2o? (auj — 


73 


= 24a^ [-1 -1- m {aoj - m)] , 


74 


= I2a^ , 


A 


= i(l 


+ 1) - 2amuj + o?uj^ + 2 . 



lim Xoc 

r* — >OQ 

um — ; — 



1 



(A5) 
(A6) 



As with the boundary conditions at the horizon, we 
construct the solution X^o beginning with boundary con- 
ditions posed at finite i?^, not at infinity. Using the 
asymptotic form of the boundary conditions to set X 
and X' at finite radius leads to errors of fractional or- 
der l/i?+ in the solution, which can be represented as 
an error in the amplitude of the out-going radiation and 
the introduction of some small in-going radiation com- 
ponent. These lead, in turn, to errors in the estimated 
power radiated to infinity and down the horizon. We 



can use Richardson extrapolation to estimate and reduce 
this error; however, as before, the radiated power must 
be determined at several different large R*^_ in order to 
estimate and reduce the error. 

To permit a more accurate estimate of the radiated 
power from X and X' evaluated at smaller we have 
solved the Sasaki-Nakamura equation analytically about 
the point at r* = oo, finding the first corrections in 1/ R\ 
to the asymptotic form of the radial function X. For 



finite R\. we have 



Xco(r*) = (l + ^) 

dXoo , ^ • /i , "^1 

— — (r*) = IH 

ar* \ r 



where 



ai 



71 
70 



uj (a^ + 4) -I- 2am -|- 



A-t-2 



(A7a) 
(A7b) 

(A7c) 



We use these expressions, evaluated at finite but large 
to set the boundary condition for the numerical so- 
lution of the homogeneous radial equation. We continue 
to use Richardson extrapolation to control the error of 
the solutions; however, each step in the extrapolation 
has a greater effect on the error and the extrapolation 
can take place at smaller R'^. 



3. A more suitable choice of variables 

The solution X to the Sasaki-Nakamura equations is a 
complex oscillatory function. Integrating the equations 
directly for X requires a spatial resolution Ar* less than 
the local wavelength of X , 



din a: 



dr* 



(A8) 



When solving for the radial function corresponding to 
a high temporal frequency |u;| the step-size can become 
quite small, with a corresponding increase in the compu- 
tational time for an accurate solution. 

It is advantageous in circumstances like these to refor- 
mulate the problem in action- angle variables, whose vari- 
ation is both slower and smoother than the variations in 
X. Writing X as 



X EE exp [i^>(r*)] 



(A9a) 



we define the two real functions ^ and <f> as the imaginary 
and real parts of <i>: 



CEE3fi($) 

= 5($) 



(A9b) 
(A9c) 



With this substitution the linear Sasaki-Nakamura equa- 
tion for complex X becomes a pair of coupled non-linear 
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equations for the real ^ and (f>. The equation for ^ is 
second order while the equation for <f> can be integrated 
immediately to obtain a first order equation. (This is ex- 
pected since the solution for X is determined only up to 
an overall phase.) Both and ^ vary slowly and smoothly 
compared to X. This is particularly true as one moves 
toward either the horizon or spatial infinity, where X is 
oscillatory in while ^ is constant and cf) is linear. Cor- 
respondingly, the numerical solution of the equations for 
4> and ^ require much less resolution for the same nu- 
merical accuracy, dramatically speeding the integration 
of the radial equation. 



4. Numerical solution of the equations for ^ and <j) 

The local errors committed by, e.g., a fourth order 
Runge-Kutta integration of the radial equation are pro- 
portional to Ar»^. Reducing the step-size and increasing 
the number of integration steps will decrease the over- 
all solution error algebraically, i.e., as a fixed power of 
Ar*, while increasing the time required for a solution. 
A higher order computational method will increase the 
solution accuracy more rapidly. Exponential convergence 
of the solution with Ar* can be obtained if the equations 
are solved via collocation pseudo-spectral techniques . 
In a collocation pseudo-spectral method the solution for 
the dependent variable is approximated as a sum over 
a suitable set of basis functions. The differential equa- 
tions, evaluated on the approximate solution at a fixed 
number of points, then determine the coefficients in the 
expansion. For problems with smooth solutions the solu- 
tion accuracy increases exponentially with the number of 
terms in the approximation (and, correspondingly, with 
the number of evaluations of the differential equation, 
which is the analog of the spatial resolution of the inte- 
gration) . Our final innovation is to solve the radial equa- 
tion using pseudo-spectral techniques. We have chosen 
a Chebyshev expansion for ^ and (j) with Gauss-Lobatto 
collocation points. Our experience is that the best per- 
formance is obtain if the integration domain RJ\ is 
divided into two parts, at approximately the peak of the 
effective potential R^: i.e., we use two expansions for (p 
and ^, one in the domain [R+,Rq\ and the other in the 
domain [Rq,R-]. At i?o we insist that the two solutions 
for (j) and ^ agree in value, and that the solutions for ^ 
agree also in their first derivative, as is appropriate for 
functions described by first order and second order dif- 
ferential equations, respectively. 
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TABLE I. n/Oisco (orbital angular velocity in units of that at the isco) as a function of r/rjaco (Boyer-Lindquist radius in 
units of that at the isco) and of a (black- hole angular momentum parameter). For a negativ e, th e hole i s coun ter-r otating 
relative to the star's orbit; for a positive it is co-rotating. This table was computed from Eqs. (3.18), (3.21) and (3.20). Near 
the isco, f2/r2isco is linear in r/nsco- 
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3455 





3480 





3536 





3572 





3664 





3885 





4069 





4611 





4930 


2 


5 





2460 





2463 





2484 





2530 





2560 





2637 





2823 





2982 





3463 





3758 


3 








1867 





1870 





1887 





1925 





1950 





2013 





2168 





2302 





2716 





2976 


4 








1210 





1212 





1224 





1250 





1268 





1312 





1423 





1520 





1827 





2025 


5 








08643 





08659 





08748 





08944 





09076 





09409 





1024 





1097 





1332 





1487 


6 








06570 





06582 





06651 





06804 





06907 





07167 





07817 





08391 





1025 





1149 


7 








05211 





05221 





05276 





05399 





05482 





05692 





06217 





06682 





08197 





09211 


8 








04264 





04271 





04318 





04419 





04488 





04661 





05097 





05483 





06745 





07595 


9 








03572 





03579 





03618 





03704 





03762 





03908 





04276 





04603 





05675 





06399 


10 








03049 





03055 





03088 





03162 





03212 





03338 





03654 





03936 





04860 





05486 



TABLE II. £ (the relativistic correction to Egw ~ ~E, the total rate of emission of energy into gravitational waves going 
both to infinity and down the hole), as a function of orbital radius r/risco and black-hole spin parameter a; cf. caption of Table 
|. This table is accurate to four significant digits; each entry was computed by summing over enough spheroidal harmonic 
orders {l,m) to produce that accuracy. 



r/n^ca I -0.99 -0.9 -0.5 fU) (h2 05 08 09 099 0.999 



1.000 


1 


240 


1 


233 


1 


197 


1 


143 


1 


114 


1 


053 





9144 





7895 





4148 


0.2022 


1.001 


1 


239 


1 


232 


1 


196 


1 


142 


1 


114 


1 


053 





9140 





7894 





4154 


0.2032 


1.002 


1 


238 


1 


231 


1 


196 


1 


141 


1 


113 


1 


052 





9137 





7894 





4160 


0.2041 


1.005 


1 


235 


1 


228 


1 


193 


1 


139 


1 


111 


1 


050 





9126 





7891 





4177 


0.2069 


1.01 


1 


231 


1 


224 


1 


189 


1 


135 


1 


107 


1 


047 





9109 





7887 





4207 


0.2116 


1.02 


1 


222 


1 


215 


1 


181 


1 


127 


1 


100 


1 


041 





9076 





7880 





4263 


0.2208 


1.05 


1 


198 


1 


192 


1 


159 


1 


108 


1 


081 


1 


025 





8988 





7867 





4434 


0.2473 


1.1 


1 


165 


1 


159 


1 


128 


1 


080 


1 


055 


1 


002 





8876 





7859 





4701 


0.2881 


1.2 


1 


115 


1 


110 


1 


082 


1 


039 


1 


017 





9706 





8726 





7882 





5182 


0.3581 


1.3 


1 


081 


1 


075 


1 


051 


1 


012 





9913 





9493 





8638 





7920 





5587 


0.4160 


1.4 


1 


055 


1 


051 


1 


028 





9919 





9733 





9348 





8583 





7960 





5930 


0.4648 


1.7 


1 


Oil 


1 


007 





9888 





9591 





9435 





9119 





8524 





8075 





6665 


0.5723 


2.0 





9893 





9862 





9705 





9448 





9312 





9034 





8530 





8171 





7117 


0.6411 


2.5 





9734 





9709 





9580 





9363 





9248 





9012 





8589 





8302 





7556 


0.7089 


3.0 





9674 





9653 





9542 





9352 





9250 





9040 





8662 





8415 





7813 


0.7469 


4.0 





9651 





9634 





9546 





9391 





9306 





9129 





8807 





8597 





8121 


0.7882 


5.0 





9665 





9651 





9577 





9448 





9371 





9216 





8930 





8742 





8320 


0.8118 


6.0 





9687 





9675 





9611 





9490 





9430 





9291 





9031 





8858 





8469 


0.8286 


7.0 





9709 





9699 





9641 





9533 





9480 





9354 





9116 





8955 





8589 


0.8416 


8.0 





9730 





9720 





9669 





9588 





9522 





9407 





9186 





9036 





8689 


0.8524 


9.0 





9749 





9740 





9693 





9607 





9558 





9452 





9246 





9105 





8774 


0.8616 


10.0 





9765 





9757 





9714 





9616 





9589 





9491 





9298 





9164 





8847 


0.8695 



17 



TABLE III. fool (the relativistic correction to -Eooi, the rate of emission of energy into harmonic-1 gravitational waves with 
frequency /i — fl/2n traveling to infinity) as a function of orbital radius r/risco and black- hole spin parameter a; cf. caption of 
Table ^. This table is accurate to four significant digits; each entry was computed by summing over enough spheroidal harmonic 
orders 2 < I < /max at fixed \m\ — 1 to produce that accuracy. 







-0.99 


-0.9 


-0.5 


0.0 


0.2 


0.5 


0.8 


0.9 


0.99 


0.999 


1.000 


3 


013 


2 


854 


2 


157 


1 


320 


1 


002 





5530 





1669 





06573 





002762 


1 


071E-4 


1.001 


3 


010 


2 


851 


2 


156 


1 


319 


1 


001 





5529 





1670 





06584 





002783 


1 


095E-4 


1.002 


3 


007 


2 


849 


2 


154 


1 


318 


1 


001 





5528 





1671 





06595 





002805 


1 


119E-4 


1.005 


2 


998 


2 


840 


2 


148 


1 


316 





9990 





5525 





1675 





06628 





002869 


1 


194E-4 


1.01 


2 


984 


2 


827 


2 


139 


1 


312 





9964 





5520 





1680 





06683 





002979 


1 


326E-4 


1.02 


2 


955 


2 


800 


2 


121 


1 


303 





9915 





5510 





1692 





06793 





003204 


1 


619E-4 


1.05 


2 


876 


2 


727 


2 


071 


1 


280 





9779 





5487 





1729 





07126 





003934 


2 


747E-4 


1.1 


2 


760 


2 


619 


1 


999 


1 


248 





9588 





5462 





1790 





07688 





005322 


5 


565E-4 


1.2 


2 


575 


2 


448 


1 


885 


1 


198 





9305 





5448 





1916 





08821 





008679 





001525 


1.3 


2 


434 


2 


316 


1 


798 


1 


161 





9111 





5465 





2040 





09951 





01272 





003073 


1.4 


2 


321 


2 


213 


1 


730 


1 


133 





8973 





5498 





2162 





1107 





01733 





005180 


1.7 


2 


085 


1 


994 


1 


588 


1 


077 





8719 





5623 





2506 





1421 





03343 





01435 


2.0 


1 


940 


1 


860 


1 


503 


1 


049 





8631 





5786 





2799 





1715 





05138 





02648 


2.5 


1 


787 


1 


719 


1 


416 


1 


023 





8599 





6045 





3237 





2151 





08223 





04991 


3.0 


1 


689 


1 


629 


1 


361 


1 


010 





8621 





6272 





3606 





2527 





1122 





07456 


4.0 


1 


567 


1 


518 


1 


295 





9987 





8704 





6638 





4191 





3143 





1664 





1222 


5.0 


1 


493 


1 


450 


1 


255 





9921 





8789 





6918 





4638 





3627 





2128 





1650 


6.0 


1 


441 


1 


403 


1 


228 





9923 





8865 





7139 





4994 





4020 





2525 





2028 


7.0 


1 


403 


1 


368 


1 


208 





9865 





8930 





7319 





5287 





4346 





2869 





2362 


8.0 


1 


373 


1 


341 


1 


193 





9829 





8987 





7469 





5533 





4624 





3168 





2657 


9.0 


1 


349 


1 


319 


1 


180 





9887 





9035 





7596 





5743 





4863 





3433 





2922 


10.0 


1 


329 


1 


301 


1 


170 


1 


005 





9078 





7706 





5925 





5072 





3669 





3159 



TABLE IV. £oo2 (the relativistic correction to Eoa2, the rate of emission of energy into harmonic-2 gravitational waves with 
frequency = 2n/27r traveling to infinity) as a function of orbital radius r/risco and black-hole spin parameter a; cf. caption of 
Table ^. This table is accurate to four significant digits; each entry was computed by summing over enough spheroidal harmonic 
orders 2 < I < Zmax at fixed |mj = 2 to produce that accuracy. 



r/rj^ca I -0.99 -0.9 -0.5 00 02 05 08 09 099 0.999 



1.000 


1 


029 


1 


020 


0.9734 





8957 





8535 





7653 





5914 





4617 





1656 





06128 


1.001 


1 


028 


1 


019 


0.9730 





8954 





8533 





7652 





5915 





4620 





1661 





06170 


1.002 


1 


028 


1 


019 


0.9726 





8950 





8530 





7650 





5916 





4624 





1666 





06212 


1.005 


1 


026 


1 


017 


0.9713 





8940 





8522 





7645 





5919 





4633 





1681 





06338 


1.01 


1 


024 


1 


015 


0.9693 





8925 





8508 





7638 





5925 





4649 





1707 





0655 


1.02 


1 


019 


1 


Oil 


0.9654 





8894 





8483 





7623 





5937 





4680 





1758 





06975 


1.05 


1 


007 





9985 


0.9548 





8813 





8415 





7587 





5974 





4773 





1909 





08241 


1.1 





9900 





9818 


0.9403 





8704 





8327 





7545 





6037 





4918 





2154 





1038 


1.2 





9648 





9574 


0.9196 





8558 





8214 





7506 





6165 





5178 





2618 





1465 


1.3 





9480 





9411 


0.9063 





8474 





8156 





7504 





6287 





5402 





3039 





1881 


1.4 





9364 





9301 


0.8977 





8427 





8131 





7523 





6400 





5597 





3417 





2277 


1.7 





9191 





9138 


0.8867 





8403 





8151 





7635 





6697 





6050 





4309 





3303 


2.0 





9138 





9092 


0.8857 





8450 





8227 





7769 





6941 





6382 





4930 





4077 


2.5 





9144 





9106 


0.8910 





8566 





8377 





7983 





7268 





6791 





5617 





4954 


3.0 





9187 





9154 


0.8984 





8684 





8517 





8167 





7525 





7097 





6075 





5524 


4.0 





9286 





9260 


0.9125 





8882 





8745 





8453 





7908 





7542 





6680 





6241 


5.0 





9372 





9351 


0.9237 





9034 





8914 





8662 





8183 





7857 





7085 





6699 


6.0 





9442 





9424 


0.9326 





9142 





9043 





8821 





8391 





8095 





7387 





7033 


7.0 





9499 





9483 


0.9396 





9232 





9145 





8945 





8554 





8282 





7625 





7295 


8.0 





9545 





9531 


0.9453 





9323 





9226 





9044 





8686 





8434 





7819 





7507 


9.0 





9584 





9571 


0.9500 





9369 





9293 





9126 





8795 





8560 





7981 





7685 


10.0 





9616 





9604 


0.9540 





9400 





9349 





9195 





8887 





8666 





8119 





7837 



18 



TABLE V. SooS, (the relativistic correction to -EooS, the rate of emission of energy into harmonic-3 gravitational waves with 
frequency f-j, = 3S7/27r travehng to infinity) as a function of orbital radius r/risco and black-hole spin parameter a; cf. caption of 
Table ^. This table is accurate to four significant digits; each entry was computed by summing over enough spheroidal harmonic 
orders 2 <l < Imax at fixed |m| = 3 to produce that accuracy. 



r/riBco I -0.99 -0.9 -0.5 00 02 05 OS 09 099 0.999 



1.000 


0.9753 





9614 


0.8926 





7848 





7309 


0.6292 





4684 





3712 





1573 





06456 


1.001 


0.9748 





9608 


0.8922 





7845 





7307 


0.6291 





4685 





3714 





1577 





06495 


1.002 


0.9742 





9603 


0.8917 





7842 





7304 


0.6289 





4685 





3715 





1581 





06534 


1.005 


0.9726 





9587 


0.8904 





7832 





7296 


0.6285 





4685 





3719 





1592 





06651 


1.01 


0.9699 





9561 


0.8882 





7817 





7284 


0.6277 





4686 





3725 





1610 





06846 


1.02 


0.9648 





9512 


0.8841 





7787 





7260 


0.6263 





4688 





3739 





1646 





07237 


1.05 


0.9507 





9376 


0.8729 





7708 





7197 


0.6229 





4699 





3780 





1751 





08391 


1.1 


0.9313 





9189 


0.8576 





7606 





7118 


0.6191 





4728 





3852 





1919 





1026 


1.2 


0.9033 





8921 


0.8365 





7476 





7026 


0.6166 





4806 





3998 





2223 





1374 


1.3 


0.8849 





8747 


0.8235 





7411 





6990 


0.6183 





4900 





4141 





2489 





1686 


1.4 


0.8728 





8633 


0.8157 





7385 





6988 


0.6223 





5000 





4278 





2722 





1963 


1.7 


0.8562 





8483 


0.8085 





7424 





7078 


0.6402 





5301 





4649 





3280 





2627 


2.0 


0.8534 





8466 


0.8117 





7531 





7221 


0.6603 





5580 





4969 





3701 





3114 


2.5 


0.8588 





8531 


0.8239 





7737 





7466 


0.6919 





5986 





5417 





4238 





3710 


3.0 


0.8676 





8627 


0.8373 





7929 





7687 


0.7190 





6324 





5785 





4657 





4156 


4.0 


0.8851 





8811 


0.8606 





8243 





8040 


0.7616 





6850 





6358 





5297 





4822 


5.0 


0.8993 





8961 


0.8788 





8477 





8301 


0.7929 





7240 





6786 





5780 





5320 


6.0 


0.9106 





9078 


0.8928 





8658 





8500 


0.8168 





7541 





7119 





6164 





5717 


7.0 


0.9198 





9172 


0.9039 





8787 





8657 


0.8356 





7780 





7387 





6478 





6045 


8.0 


0.9272 





9249 


0.9130 





8902 





8783 


0.8508 





7975 





7607 





6741 





6322 


9.0 


0.9334 





9313 


0.9204 





9004 





8887 


0.8633 





8137 





7791 





6965 





6559 


10.0 


0.9386 





9367 


0.9267 





9087 





8974 


0.8739 





8275 





7948 





7158 





6766 



TABLE VL £oo4 (the relativistic correction to Eooi, the rate of emission of energy into harmonic-4 gravitational waves with 
frequency /4 = 4i}/2Tr traveling to infinity) as a function of orbital radius r/risco and black-hole spin parameter a; cf. caption of 
Table ^. This table is accurate to four significant digits; each entry was computed by summing over enough spheroidal harmonic 
orders 2 < / < Zmax at fixed \m\ = 4 to produce that accuracy. 



r/riBco I -0.99 -0.9 -0.5 OO 02 05 OS 09 099 0.999 



1.000 


0.9393 





9209 





8319 





6981 





6342 





5196 





3574 





2720 





1116 





04673 


1.001 


0.9387 





9203 





8314 





6978 





6339 





5194 





3574 





2721 





1118 





04699 


1.002 


0.9380 





9197 





8309 





6975 





6337 





5193 





3575 





2722 





1120 





04726 


1.005 


0.9361 





9179 





8295 





6965 





6329 





5189 





3575 





2725 





1127 





04807 


1.01 


0.9330 





9149 





8271 





6949 





6317 





5183 





3576 





2730 





1139 





04940 


1.02 


0.9270 





9091 





8225 





6920 





6295 





5171 





3578 





2739 





1163 





05207 


1.05 


0.9106 





8935 





8102 





6841 





6235 





5143 





3589 





2770 





1231 





05993 


1.1 


0.8881 





8720 





7936 





6740 





6162 





5115 





3617 





2826 





1342 





07255 


1.2 


0.8560 





8416 





7709 





6618 





6084 





5109 





3697 





2947 





1547 





09588 


1.3 


0.8353 





8221 





7573 





6563 





6064 





5143 





3793 





3073 





1731 





1168 


1.4 


0.8218 





8097 





7496 





6549 





6077 





5199 





3897 





3198 





1898 





1355 


1.7 


0.8035 





7943 





7440 





6627 





6212 





5424 





4217 





3555 





2325 





1821 


2.0 


0.8021 





7934 





7493 





6774 





6398 





5669 





4523 





3880 





2676 





2190 


2.5 


0.8109 





8036 





7664 





7037 





6710 





6052 





4977 





4355 





3163 





2683 


3.0 


0.8232 





8169 





7843 





7286 





6984 





6383 





5365 





4759 





3570 





3086 


4.0 


0.8466 





8416 





8151 





7687 





7432 





6905 





5983 





5409 





4233 





3737 


5.0 


0.8656 





8614 





8389 





7987 





7765 





7297 





6450 





5910 





4758 





4257 


6.0 


0.8807 





8770 





8574 





8227 





8021 





7599 





6818 





6306 





5190 





4689 


7.0 


0.8928 





8895 





8720 





8412 





8224 





7838 





7115 





6631 





5551 





5055 


8.0 


0.9026 





8997 





8839 





8546 





8388 





8033 





7359 





6902 





5859 





5371 


9.0 


0.9108 





9081 





8938 





8673 





8523 





8195 





7563 





7130 





6125 





5647 


10.0 


0.9177 





9152 





9020 





8759 





8637 





8331 





7737 





7326 





6355 





5889 



19 



TABLE VII. Eh/Eqw (the ratio of the energy radiated down the hole to the total energy radiated) as a function of orbital 
radius r/risco and black-hole spin parameter a; cf. caption of Table |. This table is accurate to three significant digits; each 
entry was computed by summing over enough spheroidal harmonic orders {l,m) to produce that accuracy. 





-0.99 


-0.9 


-0.5 


1 


000 





0129 





0118 





00757 


1 


001 





0129 





0118 





00753 


1 


002 





0128 





0117 





00750 


1 


005 





0127 





0116 





00740 


1 


01 





0124 





0114 





00723 


1 


02 





0120 





0109 





00691 


1 


05 





0107 





00975 





00606 


1 


1 
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TABLE VIII. Af (the relativistic correction to il^ /tl = d^/dlnQ, the number of radians of orbital inspiral per unit fractional 
change of orbital angular velocity), as a function of orbital radius r/risco and black- hole spin parameter a; cf. caption of Table 
|. This table is accurate to four significant digits, and it was computed using Eq. ( |3.19| ) and £ from Table 0. Near the isco, 
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TABLE IX. T (the relativistic correction to T, the time remaining until the isco is reached) as a function of the orbital 
radius r/naco and black-hole spin parameter a; cf. caption of Table ^ Near the isco, T ~ (8/5)A/'orb oc (f — fisco)^ oc {fl — flieco)^ ■ 
We think this table is accurate to about 1 part in 500, except at r/risco ^ 1-2 where the accuracy is about 1 part in 100. The 
table was computed from Eqs. (|j.l8|), (|3.19|), and (|3.5|), using a cubic interpolation to £(r/risco) as given in Table 
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TABLE X. TVorb (the relativistic correction to A'^o^b, the number of orbits remaining until the isco is reached) 
as a function of orbital radius r/risco and black-hole spin parameter a; cf. caption of Tabic ^ Near the isco, 
T ~ (8/5)A/'orb oc (f — fiaco)^ OC (57 — fiisco )^ • We think this table is accura te to ab out 1 part i n 50 0, except at r/risco ^ 1-2 where 
the accuracy is about 1 part in 100. The table was computed from Eqs. (3.18), ( [j.l9| ), and ( |3.6| ), using a cubic interpolation to 
£{r/risco) as given in Table 
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